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SYNOPSIS 


fh© concept of ’ferkov process is used to derive 
expressions for an eduirslaat of two units in series and 
parallel considering two state fluctuating nonaal arui 
Btoray weathers, fhese expressions are then used to 
simplify a system bj series and parallel block reouction 
techaicae. It is observed that this technidue alone ie 
not sufficient for a system where stane oomplex configura- 
tion is involved. For buch cases Star-Delta conversion 
in addition with the series end parallcQ. equivalents is 
used, fhe results obtained by this technique are compared 
with those obtained by existing melhods. Finally a 
computer flow chart for system reduction is described. 
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- I 


IHTRQri/atIQK 


?erii&ps one of the most ohallerjgiag problems in this 
deoade is the design ©ad development of large scale ©yat^ss 
for coffirsercial purposes. Power sjfsteaj is also one of those 
l®rg© scale ctaasercisl syateos. In this ©g© ©nd in the ages 
coming the haaan race will be more and more depciidcnt upon 
electricity/', fhie fact introduo©© the need for reliable 
Supply to oonsumora. Historically reliability has always 
been considered during system design. However, as systems 
have become increesinisly complex the reliability problem 
has becc« 3 ® more acute. Sluoh of the early reliability work 
wm confined to making trad© offs between certain perfortaacee 
and reliability aspects of electronic aooponaats. i**or ©xaraple, 
aa the need for high power tronsnitters aeveloped, it was 
recognised that increased perfomance could only be attained 
at the price of decreased reliability. Obviously increase In 
power means greeter complexity aixi thus recuoing the mean 
time to of the transmitters. 

Uhere are various definitions of system reliability 
depending upon the particular aystsDS in which they ©re going 
to be used. However, for power systems one vihich is widely 
used is *'1?he aegree of assurance needed to deliver electrical 
power to the load, in other words availability of electrical 
power** . 



Basically there are two way® to achieve the reliability. 
The first is to develop highly reliable part® for use In eofuip- 
isent^ and The second is to design reliable systira® 

from less reliable parts. The second approach looks to be suite 
appropriate because it loads to the iacreaae in syateai reliabi- 
lity as a whole. 

The scope of this study is to develop a new technique 
for the evaluation of reliability for a complex power syoteia. 

|sl though reliability aspects sre not new their application to 
power aystsDS has started only recently • . It ia found tiiat 
the existing me'Uioda are either too compliant so or lees accurate 
for the evaluation of reliability measure®. Therefor© this 
study an effort to derive a procedure which is reasonably 

accurate and am be easily applied to even complex systems. 



- II 


^ jv* ' ?4h. *-^4 A 4n ' X it’' 


2*1 Tstyoiii^sttOi^ 

Til® reliability ftmatl<as ia tlefiaed bs tiae probebili"^ 
of the e¥«!iit naeeess for ® epeeifiei of tia«* It it 

obtaiiiesl lay mms^atlon or iat^r&tioii over that porticat of the 
e»ple spuee ^hioh defioee the oooareooe of euoceas* fMe 
iltftaition of reliability elsoes thet the relifahllitF fuwstlan 
ie directly r«lot«a to the probability fiimtiCHa* Hsuee a 
QQ^ of probabilt^ is needed to sth^ the rellebi- 

iity of a 

There ere t*o c^cz-tiZ: t'.re asessoree of prob^ility. 

First is s-priori probability or elsasical probability* Her© 
all the pooelble outcomes of m ©rperiaeiit are know* For 
crssple# whes m ooia i» tossed it ia already kncwa that it 
®ill be either head or tail sM probability of ©soh of these 
outoofti© is C*3* The etoocd is the ©xperiaoabal probabHity* 

Here the prediotions ©boat the outao^e ere based am pmt 
erperleooe* For ©xaeplej^ if sOBctMog hm imppmm^. n tiaes omt 
of S trialsf the probability of -fee eaooese f'*' TMs predto- 
tioa bec«wies better and better m H t©»d® to iafiiadty* It is 
only this R^^peot of probability efcioh is used for syetesi 
reliability etiuii*^* 

2.2 :i ecd tisT Ftul/ of .cr Reliability 

Au iooreasii^ esouot of att^atioo ia being fooue.ed on the 
i'ellability or ccaatimiity of o'ur»j.ce afforded by tmaatissloa. 




distribation and genaretion systic®. ?ise inoreesad attisntiois 
primarily stems from (i) the need to supply improved aervioe 
88 cuotc8a«*® Peowe more md more depeiaSent on their eleotrie 
service and (ii) the desire to use nm oystom voltages ai^ 
design® whose reliability is not well known to supply heavier 
loads of the future. Seeing the need of study of reliability 
one o«n make reliability models of big power ays terns end then 
try to msximize the reliability by putting more redundancy in 
the system and at the same time keeping the cost a minimum. 

2,3 Sonaeat of Power System Heliabllilv 

From a 8/steis reliability viewpoint the problem of 
prediction requires the following stepsi 

(1) Building a mfc'feemeticcl model snd defining the reliability 
measure®. 

(2) Bevsloylng a technique for ©vi^lustion of above reliability 
measures. 

(3) Comparing the predicted date and the exp erin eat al data. 

fo do this first, decision rule has to be established and 
©ecosdly effectiveness of each in terms of some physical attri- 
bute has to be exoresaed. The procedure involves developing 
acthematioel model® which sinwlate the possible outcome of 
each eltemfit® design policy. In this way it is posalble to 
compare and then to decide which ia the optioiua policy, 

Ifeasrures of effectiveness arc generally choam to 
repreomt well defined physical attributes of a syate© since 



they ©re asenable to aa theme tical model oailding, menlpulatioii 
and iieesurement. It has be®a established that the mein relia- 
bility measures in a aainteined system suoh as power ©ret 

(1) *.veil ability s fhe avellabiiity of a system is further sub- 
divided into following categories. 

(a) lastantaneoua Availability ; The probebility that the 
system will be available at time t. 

(b) Average Op Tinm t ’She proportion of time in e sped fied 
interval (C, T) that the system is available for use. 

(c) S teady State Availebility t The probability that the 
systijffl is available ©t time t >» cn. 

(2) Wem Time To Failure {*/TTTF) : The expected timo between 
system failure a. 

(3) rmration of Single Down Timeg t Some times it may be prcfersble 
to have higher systesi failure rate with shorter duration of 
sin^e down time incidents than a lower system failure rate 
vdth longer durotion of single cown time incidents. 

One should not think at this mommU that these are the 
only effective measures of power aystea reliability. In fact 
the measures are changeable eocordin^S to requireaents. For 
exscaple, a situation m®y arise where the cost of aystimo oompo- 
nents is very high. In such oases it is meaningful to cosbine 
the system perfonsanoe and above reliability measures into on# 
measure, i.e., the rstio of increase in component availability 
and increase in oompexamt oost above minlimns required. fitter 



this the optiialastion procedure aaay be to eohiov© a given level 
of ays tea reliability at least cost. 

Having obtained the reliability aobel and meesurea of a 
ayataost this study i® ooncerned with the second of the above 
aentioaed steps i.e,, developing e suitable technique for re- 
liability evcluation. The third step is retiier controversial 
ores ana not dealt in this thesis boccuse experimental data of 
a power ayatea take year® for their oolleotion. 

2.4 Type of Data lleouired 

The various data reoulred for study of power syotws 
reliobiilty can be, in gaoeral, oategorieed into following 
groups: 

Ct) Hor»ol and stormy weather component failure rates.* 

(2) Kormai and atoray weather component repair rates. 

(3) rar&tions of normal end stormy weathers. 

After obtaining these data a probability diatributton 
curve be plotted for eech group of data for a component for 
normal and stormy weathers. It has been found in reference 1 


* for any componont the failure rate I is very high in the 
beginning end in the end of it® life. In between this 
there ia a period when rote of f failure is nuite low and i« 
almost constant. If t-f is the duratiem of b ginning time 
and tg i» the duration of end time wh«in I; is very high and 
5 is the Intermediate time whcai P is fairly conetsnt, then 
for e practice! component 

t^. tg >> T 

The intermediate time T is only of interest and the component 
is used only in this period. 



that ell %hm above probability distributions are approjdlfflately 
expoaontial. fberefore if is probability of failure, 

i® tbi probtsbiiity of repair end is liie probability 

of a weather duration tiase t, then 

?^(t) « 1 - 

Py(t) w 1 - e*"^^ 

U s-iiod y in the above expressions are called svertge or 
expootea failure rate and repair rate of the cojsponent, x is 
oulled the expected duration of the weather. 

2.5 JIaaumptlona 

In order to simplify the reliability model of a system 
the followlt5g assumption© are made: 

(1) Hepair End weather duraticm distributions are exponential , 

(2) fhu tWf of the coroponeat is exponentially dtstributad. 

(3) Hepair rates of the components are much higher then their 
failure retes. 

(4) Ctoly one component fails at a time. 

(5) The partliel units in the ©yetem are completely redundant. 

(6) The system operates in two state fluctuating normal end 
stora^^ weathers. 

(7) Duration of aomsai weather is ouch higher than that of 
0 tor;iiy weisther. 

(8) Repair is possible in normal weather tmly. 



file utmtmd <|aito large izi 

finsl^art ere ealM m prsctieel power og'stee. It will be 
eboee in foUoeliig mmm sore aaeumptiooe are 


se4e in 'Varlooe approaobes for pre<llotim of power s]fsti« 
reXleMlit^ for furtter odt tbe probXeei* 



CHAFTKH - III 


BgSTIKQ j^PHlOAOHES 

5.1 Slroliricd i^p^rorch to *^ystgr. BellRbiltty^ 

fl3@ i^^thod disou&ued in this section is a block dii^rm 
redaction technique where series and parallel subunits are re» 
duced into approximate equlvnlents assuming m independent 
process and no repair In the stormy weaiiicr. Beosuae of these 
. 'il 3 ’. j the eveluatiem of reliability beooaseo very simple. 

(a) Series G.ystes! : Coiisiufer a systum composed of n dlscimilar 
components oonnected in series. The eeuirolent failure rate 
for ith ooaponont is given by 

^avi * H 4 S \ * l+S ^i 

The ^mbois used in tiie above expression are defined 
in the nomenolatare. This is m approximation when and 

fare very amoll compared to unity as Its usualXy the oeee in 
utility prsetioe. Th® reliability function of the ith component 
is given by 

H^{t) - * 

The systOT will operate only when all components will 
operate individually, therefore the reliability of the whole 
system will be 

HgCt) * H,(t).H 2 (t) ... i^(t) 

^ g-C^avl + ®av 2 ^ ^avn^^ 



fhe above expression gives the equivalent rat© of 


failure of the series 

n 



system 88 
®8Vi 


( 2 ) 


The average or equivalent repair rate can be obtained 
ea follows: 


Average repair time 


s 


average enual down tltae 

average au. of ffailuros'per year 


^avl ^ay2 ''' *** ^avn 


^avl 


D 


aY2 

0^ 


D 


avn 


n 


i«1 


'^avi* 


JL 


( 3 ) 


(b) Parallel Systea i The equations to be given for parallel 
aysteas are limited to two components in parallel. If there 
are more oomptaients connected in parallel they may be treated 
two at a time I i.e., two of the three parallel components are 
combined then the equivalent coobined with the third and so on. 
There are four ways in which the parallel system containing two 
oomponenta can fail. 

( 1 ) Initial failure is during normal weather and second 
failure is during normal weathers Since repair tifces are very 
short compared with tim^ between storms such that, at most, 
one weather change is likely during a repair time, syat«D 
failure rate can be derived as given in the next page. 



Failure rate 


* (loag run fraoti<»J of tise that wBfither is noraal) 
[(noxEial weather failure rate of component 1) 
(probability that a stora doe® not occur during 
repair of coasponmt 1) (probability that ctfflponent- 
2 fails during repair of ooaponmt 1 ) + (normal 
weather failure rate of coaponent 2) (probability 
that a storm does not occur during repair of 
ooapcaient 2) (probability that component 1 fella 
during repair of comocncnt 2) j 


Since 


10^ and HOg are much larger than 1 



it Can be shown 


that 

Failure rate ■* 


H 




N + S 12' 0 


1 


0 , 


<4) 


(2) Initial ft.ilure i® during nonual weather and second 
failure is during stormy weathers System failure rate as a result 
of first ooaponeat failure during normal weather end second 
componaat failure during stormy weather is 


Failure rate » (long run fraotion of time that weather is 

normal) [(noaal weather failure rate of 1) 
(probability that a storm ocours during repair of 
1) (probability that 2 fella during that storm) ■¥ 
(term giving other possibility)! 


1 


HU, 

(S,(£aEi*?i 52 ) 


I 




0 . 


0 , 



(5) 



(3) failure is during atoxay weather end second 

failure is during stonay wee then i»ssuKlng thet atom durations 
ere short coispared with the repair times, systets failure rate 
as a result of oiaponent failures durli^; the same storm is 

Failure rate « (expected noabor of stanas per year) [^(probability 
that 1 fails during a storm) (probability that 
2 fails diirlng the same storm) 4- (term giving 
other possibility) J 


K + S 
1 


jsrpCsD*) + (5Dj)(srp] 


,2 


( M- p'm> 

F + s' I *^1^2'' 


C6) 


It ccai b® noted that probability of a second ocMaponmt 
failure during a storm is same as the probability of the first 
failure because of the oeaorylesa nature of the exponential 
distribution describing storm durations. 


(4) Initial failure is during stormy weather and second 
failure is during normal weathers System failure rate for this 
cose is 

Failure rate * (expected no. of storms per year) [(probability 

1 foils during a ©toim) (probability 2 does not 
fail duripg atom) (probability 2 fails during 
repair of 1) + (term giving other possibility)] 


Sitioe 1 » SD^ 


Failure rate «* 



/Pi , 


(t) 


The equivalent failure rate for the parallel systoa 
in norael weather is the suiai:i®tlon of the posaihilitiea (1) and 
C 4 ) and , there f ore , is given by 

Similarly the equivalent ey'etem failure rote in the 
stomy weather in the sumsaation of the posaibilitiea (2) and 
(3) and is given by 

H.M Dp» 

“P - (9) 

The equivalent repair rat® for the syatesj ia given by 


Up . tl, ^ Ug 


(to) 


The equivalents for single state woather om be 
achieved by substituting D:| » « S « 0 in the above 

expressions. 

The above approach assusies independence of system 
states. This aesunptionf though makes the process vt^ry simple, 
is not valid for a practical system beosuse the ocoarance of 

on® ©vent depends upon the ocouraace of other or in other words 

■* 

the process is dependent. 3illinton et.al^ applied the concept 
of Markov process* which ccHialders the dependence of even 1b, 
for calouiation of pwer system reliability. The method is 
described in the next section. 



3.2 Apglloatioa of Markov Process to Power Syet^ Reliability 

A p3roc®88 is saM to be indepericient of the knowledge of 
the oatcoBie of eny preceding experiaent does not affect the 
prediction for the next experiment. For a Markov prooese the 
event® are dependent end the knowledge of the itmedlate past 
influences the prediction of following event. The theoretical 

development of the Markov approach is discussed in considereble 

4 5 

detail by Pelier and by Sandler partlouxarly with regard to 
processes that r,rc discrete in space and continuous in time. 
This metnoa clao considers a system which is continuous in 
time but aiscrete in speoe as is usually the oaae in power 
systems. The Markov technique can be applied to a single unit 
system as follows* 

The state space diagram for a single unit case is shown 
in Fig.1. The symbols used ere defined in the nomenclature. 

Considering exponential failure and repair distributions 
and negligible inoremmotal time interval dt, 

?^(t 4- dt) « P^CtXl-B dt)(1-n dt) + ?,(t)ia dtCl-"?* dt) 

+ F 2 (t)V dt(t-n dt) 

P^(t + dt) « PQ(t)n dt{1-D dt) + P^(t){1-o dt)(1-r» dt) 

+ PjCt) tJ»dtC1-.ia dt) 

fgCt + dt) « pQ(t)I3 dt(l-n dt) 4 P 2 (t)( 1 -ir dt)(1~n dt) 
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+ dt) 


» ?^(t) D* dt)+ PgCt) n dt) 

4- P3(t)(l-® dt)(1-0» dt) 

legiectiag higher order terms of dt aaad dividing by dt 
the abo m expreauions reduce to 

?o(t + dt) - Prt(t) 

^ ^ -(33 + n)pQ(t) + D!P^(t) + OPgCt) 

P-(t + dt) - P,(t) 

-J 2 np^(t) ^ im 4- 33*)?^(t) + O’PjCt) 

P«Ct + dt) - P«(t) 

^ S DPo(t) - (U -► n> + si?j(t) 

?a(t + dt) - ?,(t) 

-2 2 D»F,(t) ♦ n?2(t) - (U*4 ra)?3(t) 


faking the limit bb dt tending to zero in the above 
eaastioaa the differential equation in the natrlx form een be 
written aa 




-{3>i>n) 

m 

tJ 

0 


1 



n 

-(m+B* ) 

0 

0* 

X 

P-jCt) 

?i(t) 

IK 

B 

0 

-(tJ+n) 

n 


PgCt) 

1 


0 

B* 

n 

-(U’HKai) 


p^Ct) 


( 11 ) 


The etoohaatio transitional probability matrix will be 

0 12 3 


0 

1-B-n 

n 

B 

0 

p » 1 

a 

1-ia-B 

0 

B* 

2 

0 

0 

1-0-n 

n 

3 

0 

0* 

m 

l-ll’-m 


( 12 ) 



fo obtain the Esean time to failure it is not important 
whioh state the system feiloo in, only that it hid fail. Both 
states 2 sad 3 o®o be considered ea ebsorbiag states which ere 
states for systau failure, fhe truncated matrix X is obtain^ 
by deleting states 2 and 3 from the matrix P end is given by 

0 1 

X m 0 

1 


1-I>-n 

m 


n 


(13) 


If I is the identity metrlx then 


fi - xl 


-1 


1 




0 

iB-fd* 

a 


1 

n 

D4-n 


0 

1 


(H) 


Billinton hes shown that starting from st«te 0 the 
equivalent rat® of failure is 



SB* » M B*n 
ffl 4- E* + n 


(15) 


and starting from state 1 the ©Quivslent rate of failure is 



BD» 4> Bb 4- D*n 
m ♦ B 4- n 


( 16 ) 


The process of obtaining ecuivalent failure rate for e 
multiple unit systma is exactly similar. 5?ne differential 
equation for such a system will be of the form 

?*(t) « f P(t) 


(17) 



Xa the equation (t?) ?(t) is the transient probability 
vector and T is defined a® the system matrix of the system. 

The order of the above system matrix for a n unit systm is 
(2 2“ ) which becomes very large as n inore&ses end hence 

it beooaes too cussberaome to otot^iin its inverse in order to get 
the moan time to failure of the systes. Therefore tMs method, 
though accurate, is very difficult when applied to a system 
raving more number of units. 

The results obtoined by the method deocribea in this 
section and by simplified methou of ’’atton are oaloulated on a 
digital computer for typicul series and parnllel configaretlons 
and given in Table 1 . Because of aimplif!:.''ing asaumptions made 
Patton* a method gave widely differing equivalent rat© of failures 
for th® different case® considered. This clearly shows that th® 
assusjptior. of independence is not valid in practiccl power system 
In the next chapter m accurate method for obtaining the series 
and parallel equivalents baaed mi transitional probability matrix 
obtiined from Markov teohnique, is described. The procedure 
oombinee the merits of Petton*8 acid BiIlinton*e methode and can b 
employed both for a simple end complex configuretione with equal 


ease. 
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gsi^g PROCESS 


4*1 latrodo-otioa 

Th« preceding otiapter gav€ methods which were either 
aimple bat inaccurate or too di^Xioult though accurate. A 
ooiaproisxee between these two feature®, l.e., simplification 
and accuracy, 1® the main objective of thi® method. 7he pro- 
cedure Is same as the existing cisplifisa approach except for 
the cifferenoc that this method obtains series and parallel 
e<'ulveleat3 based on ICerkov process and therefore takes into 
account the dependence of the states of the system. However, 
certain degree of Independence still prevails becesuse of Idle 
assumption that the system cm be simplified by taking series, 
parallel snd other equivaients if spy, i.e. two or three 
eluents are taken st a time to foixs an equivalent single unit 
end then the equivalent oc®bined with others to fona another 
edulvdeat and thus the whole system in the end reduces to a 
single unit. 

4.2 Series end PerelXel Sdul valent® 

PreauBiiag that repair in stormy weather is not possible 
the state apace diagram for two elements is shown in fig. 2. the 
symbols used are defined in the nomenolature. Proceeding directly 
from the state space diagram the treoisitimal probsbillty matrix t 
and system matrix T for two elements are given by the following 



two equation SI 


0 l-D^-Dg-n 


^ 3 




0 l-O^-Dg-n Pg 




t-I5!J-I}*-a DJ B* 0 


0 


l-BJ-ffi BJ 

0 1-ffi 


•... (18) 

The a/atfsa matrix T oomes out to be the transpoae or the 
matrix (P-I) anb therefore, 








-U-j-Ug-n 


BJ 


12 

0 


m 0 


0 


-B|-m 0 

0 -BJ-ia 


.... ( 19 ) 

While deriving the equivalent single ualt for two units it is 


always assiffiieci that 


two unit system 


me unit system 


( 20 ) 



Th® identity eqiistion in the previous page saeans tfast 
all the indtipendcnt pareiaeters of a one unit system should be 
equal to all the similar paranetera of the two unit system. 
The ladeppiadtat parameters in the one unit system are 

(1) s The rate of failure of the unit starting fr«»B a 

state when unit is Aorking in norsaol weather. 

(2) SI™ t The rate of failure of the unit starting fram a 

state when unit 1® working in atomy weaxher. 

(3) y s The repair rate of the unit. 


Therefore if is the avereg;© rate of failure for two unit 

system starting fro© a system when both the components are 
working In normal weather and is the average rate of failure 

for the aystem aterting from & state when both the oomponmts ere 
working in etormy weather then. 


®avT * ®ar 
®avT “ ^av 


( 21 ) 

( 22 ) 


Substituting the values of and from ei^ations (15) 
and (16) into the above two equations 


« DP' ■» Ite D*n 

ffi -4- P* 4- n 

„ S3D* -»■ to 4- D»n 

Neglecting DP* in oomperlsm to (a® 4- D’n) the above two 

equations reduce to 

n to 4 D*n 

^avT * m 4- D*4-n 


and 


D 


avT 




( 25 ) 



( 24 ) 


■n* Sb ■* S*ii 
s 4 z> + a 


tmd 


Solving ©<?aatioa3 (23) aid (24) for U and B' 
(n+n) (bh-h) 

^ (jjHl)^^) (cnha)!?^^ “ ^^eirtP (si+a) 

(aS^y^') -jm ' 


( 25 ) 

(26) 


Tae abov® B and B* giv© the stoxay and aoraal weather 
failure r.^tes for the eqaiysd;.ent case unit oaso, la. the equaticas 
( 25 ) and ( 26 ) a end a are kac^ and the only unknowns are 
and !?he veluee of and are derived for a series 

and parallel systeas in #p endix 1.1 and 2.1 (given by the eoua- 
tiOiiS i^l.lf Al« 2 f A 2.5 and A 2 . 7 ). 

^fter equating the thtM independent perimeter of the 
single unit ayatecs to the equivalent repair rate of the two unit 
syutaui the average or equivalent one unit repair rat® for the 
parallel ana series qyatems oen be derived as shown in Appendix* 
1.2 and 2.2 (equations A ^.11 and A 2 . 11 ). 

While developing the above equivalents it was assumed 
that repair in atomy weather is aero, fhis assumption Is quite 
valid in aotuaX preotioe beoeua® in all practioal power aysteias 
the repair is dcme only after the storr is over. 

Using the above eouivslent® obtcined both for series and 
parallel oonflgoratiatis the average rat© of failures for asultiple 
unit systems for various percentages of stormy weather failures 
and differs t repair du nations are cslculated on digital computer 



&-a& givm in ?abXe 1. It osn be saen frosa this table that the 
rseaits obtfedned bj this oethod ar® fairly ideatioel to those 
ototeined by Bllliaton’s samara t© rsethod. This clearly shows 
thst syotOT reliability osn be aocurwtely aeasiired by taking 
inaiviaoul series ana parallel eqaivalent blocks. This proos- 
dare is q-iite efficiwQt for a muitlimit system. 

4.5 Technique for Complex G-Y^esa 

The method described in the last section cm be very 
well applied to a simple system^ but a situation mey arise 
where these eimple equivalents can not be dlrt.otly applied, A 
very practical eztanple will be ledder type of network as ehown 
in fig,3(a}« In this configuration the elements t end 2, 4 and 
5 sad so on ere neither in parallel nor in series. Because of 
the usebility of the above method to reduce such ayatesrs a new 
concept of Star-Delta convercion is developed. Such multiuntt 
systems can be reduced into one ebuivAeat unit by using Ster- 
Beltay eeries and parallel etoivalmts. The method is explained 
by toking the exmple of the eye ten shown in Fig. 3(a). In this 
systfss five unite ere connected in a ladder type of emfigura- 
tlon. ?lode8 a md 4 are input and output nodes. Since their 
position is not going to be chmged in the reduction process 
they are called fixed nodes. It osn be easily observed that 
nodes a* b and c are connected in delta md hence they can be 
converted into an equivalent star as ahown in Fig. 3(b) Using 
series equivslents the system of Fig. 3(b) can be reduced to 







1 


R«(*«chbr\ &jifck<in into 

a^ aiVwitfAJW" 6"^ SnAkm , 









the 0 ^’stiw of Pig.3(o). ?h@a eleaemte <2 3*) and {§ 4’) are 
in perallel m& henoe can he reduced as ^omi in fig* 3(d). 

Finally cmeapt of scries eQutvalents is once agiiin applied 
to get the ecuivalaaat single uinit ssyatem ahow In ?lg,3(s). 

In this «ay any other large power syst m asy b# token and reduc- 
ea to a saall single unit block* 

The expressions for star ©cuivnlents of a delta nod el in 
& i^ifigle 0t«te weather and as well as in two state weather are 
dertved la appendix: 5- The delta edulirslents of a star isocel 
can be derived uetiag the above expressions* For e two stftte 
wefther model the eaaiaptloa of Indf ptnJ o-'ioc of events has been 
t^ken into account* Otherwise the problem bee ernes quite complex 
©nd it is prmtioiLlly iapoealble to obtain cnalytlcal expresaloa 
for star equivalent of the deXts model* This ooticept of Star- 
Balts aanversitmy although quite old for electrical networks, is 
suooesafully appli<^ to & reliability aodel for the first time. 
The success of this epprojjoh can be observed in febl# 3 where 
the avert ge failure rates for the complex system of Fig*5Cs) 
are calculated using the above technique and Billiiitt»*e teohni** 
que and then compared. The results obtsined by the above two 
methods cms out to be almost Idc^ntiael* 



fABI2 2 Single Stet© Weathers Compari©<m of the accurate 
aethog and block diegrass reduction asethod using 
»tar-4elta comrersion (rat® of repair of a aingl# 
unit m 1000 repairs/ye er ) * 


Hat© of failure I EQuiveleiit rat© oft Bqui valent ret© of 

of a single | failure for th® 1 failure for the 

unit I OOTplex systan | cowplex ayaten 

(feilurea/vearli I Accurate Method i Block Blagreia Heduc^ 

I (gSluroe/year) | ti^ ^:ethod(FEii./^ 


0.1 

3.995x10"^ 

4.00x10*’' 

0,2 

1.699x10"^ 

1.60x10*' 

0.3 

3. 6^x10*^ 

3.69x10**' 

0.4 

6.399x10“^ 

6.40x10* 

0.5 

O *9998x10*^ 

1.00x10* 




TlJau! 3 fwo State Weathers Coapariaao of the accurate 


method and block diagram reduo ti on method using 
star<>d€lta oonveraicoa (S«1.3 hre., 1»2CX} hre.> 
S«tOOO rep*/yr.t I5gj^«0.5 fail./yr.). 


Percent age 
of atcoM^ 
weather 
failure 1 

i 

Bqui valent rate 

of failure for the amplirn 
ay stem ( fail ./yr. ) 

[ Aocorete Method 

1 

! Mook diagram Heductlcm 

I Method 

0.0 

t .0x10“^ 

1.0x10*"^ 

20 

1 .92x10"^ 

1.88x10"*^ 

40 

4.68x10*^ 

4.52x10'“^ 

60 

9.26x10“^ 

8.85x10*^ 

80 

1 .56x10"^ 

1.51x10*^ 

100 

2.38x10**^ 

2.24x10*^ 


Qm?Tm - ? 


At^?Lig»gIC!j TO ?H£ 
HSLi * BiLify : :; * by3I s 


5.1 Introduotloa 

file teohalciuee dl8ous8e<l in Chapters III and IV for 
eimpl© sYotems end for coiapiex systems om be solved on a 
digital octeiputer. In this chapter a brief review of the exis- 
ting raethoie la discasaed for the evaluatlcai of rail ability 
measures on computer, js modified flow chart is glvoa for the 
laetnod described in Chapter IV which can be used to evcluete 
reliobillty measure of evm oomplex multiple unit systems. 

5.2 Exlstlnjz Methods 

(a) Block Piaaram leduotlon feohnicue uaing Only Series 
and Parallel EouivalentS i lor a simple system where 

series and parallel subunits sre well defined the appllcnticm 

of oomputer Is fairly easy because the problem is only to ijadl- 

oate the series and parallel blocks. "BaX in case of complex 

systems where these blocks are not well observable the px^blem 

6 

is quite difficult, lor suoh ^sterns Esser et.al have used the 
concept of tie set end cut set to reduce a complex system into 
a simple syst®. fhia method involves following steps: 

(1) Pottaation of minimel tie sets of the complex systems. 

(2) finding the dual of the above set by interchanging the 
logical symbols . and in order to give the minimal cut 
set for the system. 



(3) Heproaeating the si^stOT into seriea bM psr®ll.®l block® 
alone with the help of the sbowe dual set. 

The nethod cm be aede clear by taking the exjiaple of 
the oaspiex syotera shown in Fig. 3(a)* In the cooplex ayst^ 
represent subunit 1 by subunit 2 by B, subunit 3 by C, 
subunit 4 by S and subunit 5 by S. IXiring the first step tie 
seta ere detcramed by fometim of the boolean expression 

f m 4.^ + A.a.E ^ r.E -I* J5.2.B (27) 

'Jhe dual of the above set will be 

T«» (A 4 B).(A + C + B),(B E),(Jk’a*3) 

which simplifies to 

A.B 4 A.O.E 4 B.l.D 4 3,E « Cut Sets (28) 

tJIhe above expression h®a a configuration as shown in 
fig, 4 and does not have c block other then a series or parallel 
block. Thus the conplex ^stsm finally reduces to a simple 
system and hmee oen be very wH solved on the computer using 
series end parallel eouivalmts only. 

(b) nelicbilitv *^.laatlon Uoirn: Markov Process : The 
computer program enn be made exactly on the lines of the method 
disoassed in section 5.2. The systemetio determination of the 
absorbing or falXed stivtes ciin be done bj> fonaaticn of the out 
sets of the syotem. A detailed discussion of this method is 
given in Hefer«3ce 7. 








5.5 Cosjoater froKTm Vlulm. Series. Parallel oiad Star-l^elta 

This is also 8 blo^k disgrsm rebaotLon oethod whose 
objective is to develop & prooeduie which can form the basis 
of efficient digital orasputer programs. The computer flow 
chart of this method is shown in fig. 5. Input date to the 
ccfflpater prograa is the topology of the syatme. By known 
topology erne means liie ^stem conflguratiaa in which all the 
nodes are well definedd end the reliability paroaeters defined 
in section 4.2 of the component® are given in the for® of 
matrices. In otner words one matrix glviig all the failure 
rotes in stormy weather, one matrix gtviiig all the failure rates 
in nonasl weather and on® matrix giving all the repelr rates of 
the componects connected betwe^o m^utm nodes should be given 
as the input data, the result of this ooraputer program will 
also be in the fozm of three laatrloes which will have all the 
elements except those between the input and output nodes. 

These will give ecui valent stormy weather failure, nonasl weather 
failure end repair rot® of the system. This computer pr^^graa, 
though more time oonanaing then the existing metihod using tie sets, 
is sore advant&geoub in respect that it proceeds directly from the 
topology of the system and hence saves the manual labour in forma- 
tion of minlffit'l tie and out sets. However, it is not more accurate 
because the JJtar-Dclta con version used in this program involves 
certain degree of approxioatim as disouseed in section 4.2. 
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COHCLUSIOH 


Mk acoaratc mthoa for the evaluatlcsi of po-sser a^f&tesa 
rellnbiiity mat emoider th© dependfcaoj of the troKaiti jnal 
probabilities. 1?h® existing ainple approach considers th& 
s^stenj to contain independent states hence is considered to be 
inyooarate. The existing accurate method is oonsiderably 
difficult beci-iuse the s^atm matrix to be used for evaluation 
of power a^ctm reliability is very large and it ie very much 
time oonsusijng to find its inverse. The (level opr^eat cxf edui- 
VLlmt units based on the riarkov f’rooess sots as a bridge bet- 
ween aojuriioy and siaplioity and thus sioplifies the accurate 
aetlxjd by introducing some additional approxiaetim. These 
©pproxinations simplify the reliability evrJLuction to s great 
extent and at the same time do not introduce much error in the 
results. The ooxicept of Star-3)@lta oonversicn used hare is 
Quite new to the field of reliability studies. It is shown to 
be very useful in the reduction of o<»pl 0 u systeras. m the 
end computer ©pplioj-tian has been discussed end a new eoaputer 
progruB for tho cvtluaticn of reliability is developed. The 
cesaputer program is <iuit® simple and uses only the topology of 
the systas as its input data. 

The problem of systms planning may be expected to beooae 
nor® and moze iinpoirtent in the future &m to require the deveiop- 
ment of mm advricoed teohnirues, A fatui« wonk in this area cfin 
b© to predict the number of mininuia ©tend by units required in a 



po'Aer plaat keeping the ooet to be sme of the inportsnt 
factors in the overall optiralsation. 





Unit C«39e 

2)» H normal weather f.>iiar© and repair rates, 

r* , U* stormy weather fellur© and repair rates. 

S « ^ ex!)eot 0 <i daration of stormy weather. 

jf * X expected duration of normal weather 

PpCts) probability of system being in rth Ltat& at time t 

(r can vary from 0 to 3). 

More then One Unit Case 

D^# l) 2 » 2>^» Bq, normal weather failure rates of compoacnts- 
^ 30 * ^AQ ^ 


» ^ 2 * * ^ 3 * 

®AO 

0,. 72, 7^, I'j, Ug, 
®AB’ ®A0 


Stormy weather foilure rates of componmts- 
1, 2, A, B, C, AB, B0 and AC. 

normal weather repair rates for ccasponents- 
1, 2, A, B, C, AB, BC and AC. 



Equivalent rate of failure. 



AP1>'?rJBIX - I 


pqa Tro rnr-:,! 


1.1 2~Statg Weataer 

(a) The MTTF of the ^stas Is obtained by truncating the 
stochastic transitional probability nsotri* for two units by 
elemuiiatlng the absorbing states I, 2, 3t 5, 6 and 7. 

In ieferenc© 3 it has been derived that IITTT for two 
cl-£2en'‘j3 in series atartlr^ frcaa state 0 is given by 


M. 


I)|4-33^«fEt-fn 


SO * <B^+152)Cls’+r'^4a)+n(t>'+l.p 

Therefore, average rate of failure, which is inverse 
of MTT?, for © series laodel starting from state *0* is given 
by 


B 


avS2 


(BJ+BJ+as+n) 


CA1.1) 


Similarly average rate of failure starting from state *4' 
is given by 

(B^+BJ+ia)+n(B^+Bp 
(B^+Bg+ffi+n) 


^avS2 


(A1.2) 


(b) Equivalent repair rate assuming repair is possible in 

normal weather only can be derived as fol lowing i 

The steady state availability in rth state is given by^ 

lY.J 



In the expresslcm (A1.3) |tl ia the determinant of 
th® s/stem matrix in which the first row is replaced 


1,1,... 1 and l^irl dctenalasnt of matrix Y in which 

the (l4-r)th oolUEm is repleced by 1,0,0,...0. 

for series case the system is available in states 0 
end 4. The availability in the state 0 is 

. . hioi 


TVi 


where 


(41.4) 


\Yi' 


and 


^i0‘ 
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1 
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0 
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0 

Bi 

0 

n 

0 
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-U^-Bg-n Ug 
Dg -Ug— U^-n 
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0 

0 
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BJ 

0 


t 

0 

0 


0 

0 

n 

0 


0 

0 

0 

n 
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m 

0 

0 

0 
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r,» 


1 

m 

0 

0 
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-B^-m 


B^ 


1 

0 

IQ 

0 

0 

0 

-BJ-ffi 
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0 

» 

0 

0 

0 
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1 

0 

0 

Bt 

0 

0 

0 

-m 

1 

0 

0 

St 

0 

0 

0 

-a 



The value of these deterainsnts is found to be 


Y « -®yi-*'a*y2*(a+n)yj+ffiCB^4*BJ+ia+n)y^-(Ji»4ffi)y^+(DJ+ai)yg (A1.5) 
where 

♦tllUgCDIj+m) (D:|+ia)4i3£^ J 

yj “ -<I>2+®)[Sa(®t+Sp^(2>|+I>2)(33^4Bj4-iB)]lJ2 [35^(15, 4n)+inD,] 

[n(B]4D* )+B^ (l>*+I)*4ai)J ^Og^® (nUj+B^ (B*4l)j4ia)] ] 

y^ • finD<j(lJ<j4l)2)4]D^BJ(U^+B24n)] [(B54ia)I)2{l5;4l}J+ia)4iBaD* ^ 

+ [[nS^4D, (B'+BJ+ai)] { (»2+n)BJ+®B2] 

4£ijDjj]+[iaI)lj<U^4B2)4l):| PjCU^+Bg+n)] [nI3J(D^+02'*‘“) 
4nB2<I>^+BJ+ffl)] 

y4 » -O^Ug [(15^40) {Bj4i») (l)^4T)2+U^"»-02)+nB2^B»4is)4aD5 (BJ4 ib)] 

4n0^ ^ [BJ [»1)^4B2 {Bij4r*j4 ia) (Ug+a)] ♦a {%( 33 j +15 J 4 a) 

jg « -nUg [B^ (D^4ii)4aD,] [al)^4B2(Bj4a)4B^ (13^4Bj-».ffl)4.I>* (I)2+U,4n) ] 
-nOj2[B^ (Ij*3'J«)+D5 (U,™)] « {l)-,(B|4l3j4ffl) 

4B| (U^4Ii)]] -aBg [0^1)^400,] [l)^CB|4Bj4ffl)4r| {B2+0,4n>J 



4.BJIl2+atJ2)+(B^+a) <B^B;<fis33^4n]}p ] 


^Yailabllity In the stste 4 Is given by 



(# 1 . 6 ) 


1^141 Is obtained by replacing the 5 th ooluain of [T 3 
by 1 » 0 , 0 t... 0 . It oojosea out to be 

1 ^ 14 ! • -tr|U2*sn[C0^+U2+S-^+B2)(l>pffl)(D^4'ffl)aP^(rpffl) 

+nI)|(T!^+ai) J (Al. 7 ) 


Hence total steady state availability in the series 
case is 

■^10 1 + 

■*82 - “b ♦ "4 - -“ - -|Y T ^ 

The steady state unavailability for aeries case is 
given by 

^ 512 **^“ %2 (^ 1 * 9 ) 


For eapoaentially distributed MTfF and repair time the 
unavailability of a ©yetem is tie ratio of eouivalent failure 
rate eoad eduivalent repair rate of the ayetemi therefore. 


rt « ^avS 2 
'*32 * “ 1 ^ 


(A1.10) 


where Bgg is the equivalent repair rate for the series model. 



Tais gives 


or 
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\vB2 




Yl D 


avS2 




(A1.11) 


1.2 Instate Weather 

1-3tste weather is eouivalent to 2-8t8te weather with 
3};j ** D2 * h ® 0 and ffi a# 1 . 

equivalent rate of failure is obtain edi by substituting 
* I>| ** n » 0 and ra ■ 1 in expression (A1.1), i.e., the 
expression for eouivalent failure rate for two-state weather 
case. Thus the eouivalent rate of failure for sin^s state 
weather model is obtaln«i as 


®av31 - ^ 

B* . DJ » n - 0 md m - 1 in (A1.11) 

the eoulvalfflat repair rate for sin^e state weather series model 
is given b/ 

'^31 - — 


a. U, , Vj » ^ 

^31 ^ 


(AI.I 3 ) 


The expressions for and are saaie as derived 

7 

by Stanton . 



A ^ ” .TTY - II 


Ea?i?j>i.i?arT modkl fob 7wo omits im ?mAmsL 

2.1 2~State Weather Considerlry. Regatr Rates in Stomy 
Weathar to be Zero 

(a) KCjujyalent Failure Rate t Accordijqg to Billinton^ 
the MTT? is obtained by truncating the stochastic transitional 
probability iiatrix by deleting the absorbing states and then 
inverting the truncated matrix after subtracting identity 
matrix from it. 

Therefore, if the truncated matrix is *X* wildi 
z transition states, then let 




®12 

e e • 

.... 

• « • # e 

(I— X) * —A * 

®r1 

e e e 

®r2 

« « • « 
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e ♦ * * • 
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The MTTF starting from rth state is given by 

(Eyi + Ky2 ^ •••• + ^1*2^ (A2.1) 

where I ie the identity matrix and ^ 3 . 2 **“* ^ra 

the oof actors of 8^,^, aj;.2***** ®r« ^ above matrix. 


It oan be eaail;^ ahowd that 
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zz 
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Define the above die term iaesxt by -IaI^* This matrix: is 
obtained by adding one row with all Is and one colyyian with all 
aeros exoopt in the rth row, where it is 1 , in the matrix -A * 

Henoe the average rate of failure starting from 
rth state is 


T. . JL - U' 

\rr My lljy 

For the present case the detem-nant 


{A2.2) 

A I comes out 


to be 


\k\ m [nD^+D^ [nD^+DgCDj+ai)] (Dij+DJ+M)+n(D^+Dp] 

[x*i(B^+33j4SB)+nB| ] 


(A2.3) 



aad 


4*B ) -fxiBg ( D^ "^^2 

♦aa+n)] + ) (B^+is)^^^;] [(B’+BJ^m) [(DJ+ai) (B^+Bg) 

+iBB^4-a(l}|4-B2)-»-B^I)|,j-i-ii(X!^+Bj|) (0^+D2'*'’^5'*‘®^(®1**‘^1'»'®2) 1 

» • • • CA2*4) 

!Pherefore average rate of failure starting from 
state 0 is 


"av?2 * pTf^ 

Starting from state 4 


(A2.5) 

I X can b. Sifflpliflad to b« 


II 4 - y^-yj+yj-^i+y^-y's 


(A2.6) 


where 

{BJ+ffl)+nI)j]+0^B^ [(tJg+D^ ) (D^^-m)+nB| ] 

^ [(B^+Bg) (IS^4m)+nBj| [ ) (B|+m)+nB5 j ] 

+I)^(Bj4-ffi)[(B2'^»i)(2)^+»)+B2DJ I J 
y| » -a [(Bg+B, ) (Bij4-ffi34*nB^ ] {(B^+Bg ) (B^+m)4tiBj|-B^|»2^^2‘^®^ 
4aB J ^ [ (^2’*h ) )+nB^ I -Bg [®1 (® t ■*'® 1 ] [^^1 

4>B2)(B^4!B)+nt}J ] 





This average rate of failure starting froa state 4 1» 


B* 


■»■-*§;., 1, 


(A2 


t) 


•*’? \j\^ 

(*>) Enul valiant Heoair Bate s Tii® systea is unavellebla ^ 


states 5 arxJ 7. Unavailability in state 3 is given by 


» 


Yi.i 


(#2. 


0 ) 


*5 ” i 1, 

where | Y ) has been already derived in Appmdix I and ] | 

be alaplified to be 


I « -a[[(ffi+1)J)(U^D^4D^B2)+nBp^]{Cl>|4a)(I»j4Bj4B)B2 

+DjDBn]]-a \n'B\ rjU^+ [nT>*4l)^ (Dlj4Bj4m)]{aI)J 

4B2(B|+ra)]] [nl)J(B-,4U24n)+nB2(B|4Bj4a) J 


Unavailability in state 7 is given by 


Y, 




17 


CJ»2,9) 


where | | is siapliflea to be 

I '^17 1 *• 3*^1 y 6 

where end yg are already obtained in Appttidijc I. 



Hence total unaYallebillty for parallel eyst®a 
is given by 

I . I i 


+ 5, - llislliljl! 

“p "3 * ’7 trr 


avPg 


0 


P2 


where Cfpjg is the equivalent repair rate for two state weather 
nod cl. 


Trom the above expression 


D 


m 


avP2 


^ ■ ' ^i5rji7i 


(A2.10) 


2.2 Single State Weather 

Substituting * n « 0 and a « 1 in equation (A2.5) 

the equivalent rate of failure for single stete unvlroma^t 
parallel nod el is given by 

®avP1 “ "(y, •►Bg ) j 


Since B^, Bg 

B,B2(tJ,4.0g) 

^av?1 S^Ug 


(A2.11) 


fhe equivalent repair rat® can be obtained by re- 
subet ituting B| * BJ « n * 0 and ns » 1 in oquaUcm (A2.9)* Heno 

n . (n,*ll2)(0,+I)i)(Vl>2) 

(tl^+Ug) (A2.12) 

7 

These expressions are s«e as those derived by Stanton * 



APPMDIX - III 


s?AR~i)EL*gA oomrmsTQir 


CoQsider elanoats AB, B0 ©ad AO ooaatoted in delta 
liave an equivalent star with elementa A, B and 0 as shown 
in Fig* Al* The delta equivalent of the star nodal for 
single state as well as for two state weather case la 
obtained In the following sections; 


3*1 Single State Weather 
In the Fig.Al 

A and C in series s. (AB and BO in series}in parallel with AC 


Using equations (A1.12), (Al.13)t (A2.11) and (A2. 12) 

'AC 


rs ^ n fr^ \-n fJi , ^AB^BC'*‘^BC^AB 1 


(r X - :V - [®Ab''bo(“aB*®AO>* ] <«• ’ > 


AO AB 30 


Si.iilarly by syametiy the expressions for (B and C) and 
(A and B) in series are 

®B ♦ ®o - [®ab<'’ab*“bo>®a5*®ao“ab«’»bc*''ab)] <«-2) 

fU fin i3v 

- Ko'’bo(“ab*’»ac>*®bo"ao(''ab"’'bo>] 

AO AB 30 


\ ^ % 




Equations giving equivalent repair rates ai 


«a'’o<V®o) 

^wW 


“ab^bo^^bs^ab 


V M ,.v. . 0 „ 


« tJt— 


AO'^ATJ 


Vb<V^b> 


iy.« -♦• 

AB 


^bc^'acC^bc^^ac^ 


Solving equations (A3»1) to (A3«6) expreasimi 
equivalent star feilare and repair rates are 


», - < S >< • 


'3 


°AB°3C )( „ +„ ) 

A 5 .: V 


B, 


0 




and 


B. 

A 

”aC AB 

Ul 

"b 

« 

BO AB 

U3.i 

Bo 

* ^AC ^30 

U3.1i 


3*2 Two State Weather 

It is very difficult to get the expresslcms for star 
equivalent for a two state weather laodel using series si^ 
parallel equivalents obtained by Markov technique. In order to 



simplify the process following approjdjaations are made 


(1) ^ is very srasll (i.e. 1 » ^ ) 

(2) The process is independent 

(3) Repair is possible only in normal weather. 

The lest assumption has been already used for a two 
state weather model. 

1 

f>erit,a equivalent for two elsmente is given by 


and 


or 


D* W D! + Di 
e :5 -^2 

HCD^+Bg) + S(I)^ 4 B*) 


(A5.t3) 

(A3.U) 

(13.15) 



Since 1 f 



4 . 



(B,4l]^) 



1 




The second term is very amall* therefore, 



for two elesjtnts in parallel equlTraleat 


l>eP “ 


4. -5 

^2 

(A3.17) 

®;p ■ 

B-Bl 

u, * -§7* 


(A3.18) 

''.p - 

C, ♦ Hg 


(A3.19) 


Following the saise procedure as in App«aadix 3.1 the 
eouationa for (A and a)i (A and 3) and (B and C) in aeries 
giving the nonaal weather failure rates can be written as 




T)* D* 

+ ^\] 

'Vw ^ U "***/ J 


(A3.20) 




/|»U 4.gT)t fcif.JlS. 

AB ^ 


. ^)] 


<A3.21) 




AB ''AG 


.S(^ * ^)] 

^AB AC 


CA3.22) 



On solution of equations (A3. 20), (A3. 21) and (A3. 22) 


®A " 


_ _ / 1 . 1 %. 


'AB 


Bg 


-n f J- X JLu 4. 

AB BC tJ^g UjjQ I '^1^3 


{A3.23) 


(A3.24) 


„ _ ,1 . 1 %. S/^AC^io . %C®AOx 

j>o " °ac»do< 5^ * s'-TJiT * 

iquatlona giving storeiy weather failure® are 

B'.Bi . 


(A3. 25) 


(A3. 26) 
(AI.27) 

(A3.28) 


BJ 






Solving equrations (A3. 26) to (A3. 28) 

. ♦ ^s^ab^bo 

^AB BO 

, ^AO^fc . ?io?BC ^ n« . 


{A3.29) 

(A3.30) 

(A3.31) 


^AO BO 

Since single and two state weather equivalent repair 
rate® for aerie® and parallel model come out to he approacimstely 
same, the equations giving repair rate® for the equivalent 


star model for two state weather cose will be seme as the 
equations (A3.4)i (A3. 5) and (A3. 6). 


The equations (A3. 4), (A 5.5) and (A3.6) can be solved 
t 1 1 

for IJ- » nr ^ written in the following matrix fonsii 

A B C 


' 1 " 


1 

1 

1 

A 


57 

A 

’®A 

h 

1 

1 

1 

-L 

JL 


* S' 

2 


Bg 

^3 

1 


1 

1 

1 



1 


1 

w| 
o 1 

1 


%B^BO ^^AB‘*'^AC ^ C'^Sc ^ 


^AS^AC (^/ )^®A O^AB (^AC'*‘Sc ) 




(A3.32) 


Substituting the values of Bg and Bq for 
equations (A3. 23) i (A3. 24) and (A3. 25) and maxing the approxi- 
mation that 1 » ^ t it can be shown that 


"a 


(A3. 33) 

"b 

* ®AB “bo 

(A 3. 34) 

“o 

" “ao * “bo 

(A3. 35) 
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